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ABSTRACT 

We show how the complex-shift method introduced by Appell in gravity to the 
case of a point mass (and applied among others in electrodynamics by Newman, 
Carter, Lynden-Bell, and Kaiser to determine remarkable properties of the electro- 
magnetic field of rotating charged configurations), can be extended to obtain new and 
explicit density-potential pairs for self-gravitating systems departing significantly from 
spherical symmetry. The rotational properties of two axysimmetric baroclinic gaseous 
configurations derived with the proposed method are illustrated. 

Key words: celestial mechanics - stellar dynamics - galaxies: kinematics and dy- 
namics 



1 INTRODUCTION 

For the discussion of many astrophysical problems where 
gravity is important, a major difficulty is set by the po- 
tential theory. In general, to calculate the gravitational po- 
tential associated with a given density distribution one has 
■ to evaluate a three-dimensional integral. Except for special 
circumstances, where a solution can be found in terms of ele- 
mentary functions, one has to resort to numerical techniques 
and sophisticated tools, such as expansions in orthogonal 
functions or integral transforms. 

Under spherical symmetry, the density-potential rela- 
tion can be reduced to a one- dimensional integral, while for 
axisymmetric systems one in general is left with a (usually 
non-trivial) two-dimensional integral. As a result, the ma- 
jority of the available explicit density-potential pairs refers 
to spherical symmetry and only few axially symmetric pairs 
are known (e.g., see Binney & Tremaine 1987, hereafter BT). 
In special cases (in particular, when a density-potential pair 
can be expressed in a suitable parametric form) there exist 
systematic procedures to generate new non-trivial density- 
potential pairs (e.g., see the case of Miyamoto & Nagai 1975 
and the related Satoh 1980 disks; see also Evans & de Zeeuw 
1992; de Zeeuw & Carollo 1996). For non-axisymmetric sys- 
tems the situation is worse. One class of triaxial density 
distributions for which the potential can be expressed in a 
tractable integral form is that of the stratified homeoids, 
such as the Ferrers (1887) distributions (e.g., see Pfenniger 
1984, Lanzoni & Ciotti 2003, and references therein) and 
special cases of the family considered by de Zeeuw & Pfen- 
niger (1988). Additional explicit density-potential pairs are 
given by the Evans (1994) models and by those constructed 



with the Kutuzov-Osipkov method (1980, see also Kutuzov 
1998). Another way to construct in a systematic way explicit 
density-potential pairs with finite deviations from spherical 
symmetry has been presented in Ciotti & Bertin (2005). 
This technique is based on an elementary property of the 
asymptotic expansion of the homeoidal potential quadra- 
ture formula for small flattening^], and recently, it has been 
applied to the modeling of gaseous halos in clusters (Lee & 
Suto 2003, 2004), to the study of the dynamics of elliptical 
galaxies (Muccione & Ciotti 2003, 2004), and finally to a 
possible interpretation of the rotational field of the extra- 
planar gas in disk galaxies (Barnabe et al. 2005, 2006). 

In the context of classical electrodynamics a truly re- 
markable and seamingly unrelated result has been obtained 
by Newman (1973, see also Newman & Janis 1965, Newman 
et al. 1965), who considered the case of the electromagnetic 
field of a point charge displaced on the imaginary axis. Suc- 
cessive analysis of Carter (1968), Lynden-Bell (2000, 2002, 
2004ab, and references therein; see also Teukolsky, 1973; 
Chandrasekhar, 1976; Page, 1976), and Kaiser (2004, and 
references therein) revealed the astonishing properties of the 
resulting holomorphic "magic" field. As pointed out to us by 
the Referee, the complex-shift method was first introduced 
by Appell (1887, see also Whittaker & Watson 1950) to the 
case of the gravitational field of a point mass, and succes- 
sively used in General Relativity (e.g., see Gleiser & Pullin 



1 Such expansion can be traced back to the treatise on geodesy 
by Sir H. Jeffreys (1970, and references therein; see also Hunter 
1977). 
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1989; Letelier & Oliveira 1987, 1998; D'Afonseca, Letelier & 
Oliveira, 2005 and references therein). 

In this paper we show how the complex-shift method 
can be also applied in classical gravitation to obtain exact 
and explicit density-potential pairs deviating significantly 
from spherical simmetry; we are not aware that this possibil- 
ity has been noticed in the literature. The paper is organized 
as follows. In Sect. 2 we present the idea behind the method, 
and we discuss its application to the case of spherical sys- 
tems, describing two new self-gravitating and axisymmetric 
systems obtained from the Plummer and Isochrone spheres. 
In Sect. 3 we focus on the rotational fields of the new pairs 
when interpreted as gaseous systems. Section 4 summarizes 
the main results obtained, while in the Appendix we present 
the peculiar behavior of the singular isothermal sphere un- 
der the action of the complex shift. 



2 GENERAL CONSIDERATIONS 

We start by extending the complexification of a point charge 
Coulomb field discussed by Lynden-Bell (2004b), to the 
gravitational potential "^(x) generated by a density distri- 
bution p(x). From now on x = (x,y,z) will indicate the 
position vector, while < x, y >= Xii/i is the standard inner 
product over the reals (repeated index summation conven- 
tion implied). 

Let assume that the (nowhere negative) density distri- 
bution p(x) satisfies the Poisson equation 



A-xGp, 



(1) 



and that the associated complexified potential "1> C with shift 
ia is defined as 



$ c (x) = $(x — ia) 



(2) 



where i 2 — — 1 is the imaginary unit and a is a real vector. 
The idea behind the proposed method is based on the recog- 
nition that 1) the Poisson equation is a linear PDE, and that 
2) the complex shift is a linear coordinate transformation. 
From these two properties, and from equations (1) and (2) 
it follows that 



V J $ C = 4nGp c , 
where 

p c (x) = p(x - ia) 



(3) 



(4) 



Thus, by separating the real and imaginary parts of $ c and 
p c obtained from the shift of a known real density-potential 
pair one obtains two real density-potential pairs. 

A distinction is in order here between electrostatic and 
gravitational problems: in fact, while in the former case a 
density (charge) distribution with negative and positive re- 
gions can be (at least formally) accepted, in the gravita- 
tional case the obtained density components have physical 
meaning only if they do not change sign (see however Sect. 
4). Quite interestingly, some general result about the sign 
of the real and imaginary parts of the shifted density can 
be obtained by considering the behavior of the complexified 
self-gravitational energy and total mass. In fact, from the 
linearity of the shift, it follows that the volume integral over 



the entire space 
W c = ~ j p c $ c d 3 x 

= \ I [K(p c )K($ c ) - 9(p c )9($ c )]d 3 x 



(•») 



coincides with the self-gravitational energy W = 
0.5 J p"I > d 3 x of the real unshifted seed density. Thus the 
imaginary part of W c is zero, i.e. 



- J [K(p c )9($ c ) + 9(p c )K($ c )Kx = 

-G I / M^dVx^O, 
x — x' 



(6) 



and Wc is the difference of the gravitational energies of the 
real and the imaginary parts of the shifted density. The van- 
ishing of the double integral (0 shows that the integrand 
necessarily changes sign, i.e., the complex shift cannot gen- 
erate two physically acceptable densities. Additional infor- 
mations are provided by considering that, for the same rea- 
sons just illustrated, also the total mass of the complexified 
distribution M c = J p c d 3 x coincides with the total (real) 
mass of the seed density distribution M — J pd 3 x, so that 



3?(p c )<f x = 0. 



(7) 



Identities (©-([7} then leave open the question whether at 
least K(p c ) can be characterized by a single sign over the 
whole space. It results that either cases can happen, depend- 
ing on the specific seed density and shift vector adopted. In 
fact, in the following Section we show that simple seed den- 
sities exist so that 5ft(p c ) is positive everywhere, while in 
the Appendix we present a simple case in which also 5R(p c ) 
changes sign. 



2.1 Spherically symmetric "seed" potential 

We now restrict the previous general considerations to a 
spherically symmetric real potential $(r), where r = ||x|| = 
y/< x, x > is the spherical radius, and ||...|| is the standard 
Euclidean norm. After the complex shift the norm must be 
still interpreted over the realfl (e.g., see Lynden-Bell 2004b), 
so that 



2 n . . . 2 
r — li < x, a > —a , 



(8) 



where a = ||a|| . Without loss of generality we assume a : 
(0, 0, a), and the shifted radial coordinate becomes 



Harp — a 



p = cos f 



(9) 



where 9 is the colatitude of the considered point, pr = z, 
and 



Pc 



Haz — a 2 



2 — 2iaz — a 2 



(10) 
(11) 



If one adopt the standard inner product over the complex field, 
one would obtain llx — ia|| 2 =r 2 + ||a|| 2 . 
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Note that, when starting from a spherically symmetric seed 
system, the real and immaginary parts of the shifted density 
p c can be obtained 1) from evaluation of the Laplace opera- 
tor applied to the real and imaginary parts of the potential 
$ c , 2) by expansion of the complexified density in equa- 
tion (1111) . or finally 3) by considering that for spherically 
symmetric systems p = p(<&), and so the real and imagi- 
nary parts of p c ($ c ) can be expressed (at least in principle) 
as functions of the real and imaginary part of the shifted 
potential. 

2.2 The shifted Plummer sphere 

As a first example, in this Section we apply the complex shift 
to the Plummer (1911) sphere. We start from the relative 
potential ^ = — $, where 



GMl r- 



(12) 



and r is normalized to the model scale-lenght b, so that the 
associated density distribution is 

"=4^3? ( 13 ) 

(e.g., see BT). For ease of notation, from now on we will use 
normalized density (to M/b 3 ) and potential (to GM/b), and 
so 



(14) 



From substitution (8) (where also the shift a is expressed 
in b units), the shifted potential \& c = l/£ e depends on the 
square root of fj? = 1 — a 2 + r 2 — 2iaz = de bv> , with 

d = |Cd 2 = v/(l-a 2 + r 2 ) 2 +4a 2 2 2 , (15) 
and 

1 - a 2 + r 2 2az 



cos(p= , sin<p = — . (16) 

Note that cos cp > everywhere for a < 1, and in the follow- 
ing discussion we restrict to this case. The square root 

C c = Vde 71 ^^ 2 , 



(* = 0,1), 



(17) 



is made a single-valued function of (r, z) by cutting the 
complex plane along the negative real axis (which is never 
touched by £jf) and assuming k = 0, so that the model equa- 
torial plane is mapped into the line ip — 0. With this choice 
the principal determination of c reduces to \P when a — 0, 
and simple algebra shows that 

\/2az 



cos — = 



y T+ cos <p 

V2 ' 



• <P 

sin — = — - 



2 s/2 2 dy/l + cos ip 

The real and imaginary parts of ^ c are then given by 



Co 
ICcI 2 



Vd + 1 - a 2 + ; 



V2d 



(18) 



(19) 



(20) 



respectively, and from equation (|14|l we obtain the expres- 
sions of the (normalized) axysimmetric densities: 



47T 



3f(* c 



10a 2 z 2 
d 4 



+ 



d 8 5R(* c 



(21) 



3S(^ C ) 

47T 



55R(* C ) 4 - l^jl- + 



d 8 3?(* c ) 4 



(22) 



We verified that the two new density-potential pairs (|19l) - 
pl]l and p0l) - ((22)) satisfy the Poisson equation {TJ. 

Note that, at variance with 5R(»1/ C ), the potential £s(*l/ c ) 
change sign crossing the equatorial plane of the system. In 
accordance with this change of sign, also Q'(pc) is negative 
for z < 0, as can be seen from equation (122[) . and thus 
cannot be used to describe a gravitational system; we do 
not discuss this pair any further, while we focus on the real 
components of the shifted density-potential pair. Near the 
centerd~ 1 - a 2 + r 2 [l + 2a 2 p 2 /(l - a 2 )] + 0(p 2 r 4 ), and the 
leading terms of the asymptotic expansion of equations (|19|) 
and (1211) are 



5ft(* c 



5R(p c 



r 2 (l-a 2 +3aV) 4 

2(1 _ a2)S /2 (23) 

15, 2 (l-a 2 +7aV) +0(r4);(24) 



47r(l-a 2 ) 5 / 2 87r(l-a 2 ) 9 / 2 



in particular, the isodenses are oblate ellipsoids with minor- 
to-major squared axis ratio (1 — a 2 )/(l + 6a 2 ). For r — > oo, 
d ~ r 2 + 1 - a 2 + 2a 2 /i 2 + 0(/i 2 r _2 ), and 



»(^ e )~i- j-^+ggV +o( 

r 2r a 



(25) 



15/1 1 — a 2 + 7a 2 u 2 \ 



2r 2 



so that 5R(p c ) coincides with the unshifted seed density (13) 
and it is spherically symmetric and positivfl Thus, near 
the center and in the far field 3?(p c ) > for < a < 1. 
In addition, on the model equatorial plane z — (where 
d = 1 — a 2 + ii 2 , and 7? is the cylindrical radius), 5R(^1/ C ) 
coincides with the potential of a Plummer sphere of scale- 
lenght vl — a 2 , and from equation l|2ip it follows that 
K(p c ) = 3K(* c ) 5 /(4tt) > for < a < 1. 

However, for z/Oa negative term is present in equa- 
tion (|2ip . and the positivity of 5R(p c ) is not guaranteed for a 
generic value of the shift parameter in the range < a < 1. 
In fact, a numerical exploration reveals that 3£(p c ) becomes 
negative on the symmetry axis J? = 0atz~0.81 for 
a = a m ~ 0.588; the negative density region then expands 
around this critical point for increasing a > a m . The iso- 
density contours of 5t(p c ) in a meridional plane are shown 
in the top panels of Fig. 1, for shift parameter values a = 1/2 
(left) and a = 23/40 = 0.575 (right). The most salient 
property is the resulting toroidal shape of the model with 
the large shift, which reminds similar structures known in 
the literature, as the Lynden-Bell (1962) flattened Plum- 
mer sphere, the densities associated with the Binney (1981) 
logarithmic potential and the Evans (1994) scale- free po- 
tentials, the Toomre (1982) and Ciotti & Bertin (2005) tori 
(see also Ciotti, Bertin, & Londrillo 2004), and the exact 
MOND density-potential pairs discussed in Ciotti, Nipoti & 
Londrillo (2006). Unfortunately, we were not able to find 
the explanation (if any) behind this similarity (see however 
Sect. 4). 



3 This is a general property of shifted spherical systems, as can 
be seen just by expanding equation lllll l for r — * oo. 
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Figure 1. Isodensity contours in the (R,z) plane of 5R(p c ) of the shifted Plummer sphere for a = 1/2 (top left) and a = 23/40 (top 
right), and of the shifted Isochrone sphere for a = 1/2 (bottom left) and a = 4/5 (bottom right). The coordinates are normalized to the 
scale-lcnght b of the corresponding seed spherical model. 



2.3 The shifted Isochrone sphere 

Following the treatment of the Plummer sphere, we now 
consider the slightly more complicate case of the shifted 
Isochrone sphere. Its relative potential and density are given 
by 



P = 



GM 



b 1 + C 
M 3(l + C) + 2r 2 

4^b~z (i + C) 3 C 3 ' 



(27) 



(28) 



where £ is defined in equation (12), and again all lenghts are 
normalized to the scale b (Henon 1959, see also BT). In anal- 
ogy with equation (14), the normalized shifted density can 
be written as a function of the normalized shifted potential 



':! • 2r;M'.. 



4tt<, 



The real and imaginary parts of $ c = (1 + £ c )/[l + \Q C \* + 
25f(£c)] are easily obtained from equations (15)-(18) as 



(29) 

(l + Cc)/[l+ICc 12 



3(* c 



1 + Vd+l-a 2 +r 2 /V2 
l + d+y/2(d+l-a 2 + r 2 Y 

\/2az 



(30) 



y/d(d+l-a 2 +r 2 )[l + d+ ^2(d+l-a 2 +r 2 )\ 



while 5R(p c ) and S?(p c ) can be obtained by expansion of 
equation (29); however their expression is quite cumber- 
some, and so not reported here. Again, Ss(\l/ C ) changes sign 
when crossing the model equatorial plane, revealing that 
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also 3(p c ) changes sign in order to produce the resulting 
vertical force field near z = 0, and the vanishing of the in- 
tegral in equation ([7|. In accordance with equation (llip on 
the model equatorial plane 5R(p c ) is functionally identical to 
an Isochrone sphere of scale-lenght \/l — a 2 (and so positive 
for < a < 1), while for r — * oo 



[R(*c 



r - 1 1 



2r 3 



o 2 2 



(32) 



^ ~ 4^ ( 2r - 3 + 4a2L V^j + °(^ 7 )- < 33 ) 

As for the shifted Plummer sphere, also in this case 5ft(p c ) 
coincides in the far field with the seed density, and so it is 
positive for < a < 1, while it becomes negative on the z 
axis at z ~ 0.648 for a > a m — 0.804. In the bottom pan- 
els of Fig. 1 we show the isodensity contours in the merid- 
ional plane of 5ft(p c ) for the representative values of the shift 
parameter a = 1/2 (left) and a = 4/5 (right): again, the 
toroidal shape and the critical regions on the symmetry axis 
are apparent in the case of the larger shift. 



R 



I dp 9$ dp 9$ \ 
\dRdz 7 ~ Ite'dR.) Z 



(37) 



We remark that this "commutator-like" relation is not new 
(e.g., see Rosseland 1926, Waxmann 1978, and, in the con- 
text of stellar dynamics, Hunter 1977). 

Before solving equations (|34|l for the real parts of the 
density-potential pairs of Sect. 2, it is useful to recall some 
basic property of rotating fluid configurations. For exam- 
ple, in several astrophysical application (e.g., the set-up of 
initial conditions for hydrodynamical simulations), equa- 
tions (|34[) are solved for the density under the assump- 
tion of a barotropic pressure distribution (and neglecting 
the gas self-gravity), and for assigned p(R, oo) = or 
p(R,Q) = po(R). 

As well known, this approach can lead only to hydro- 
static (Q — 0) or cylindrical rotation (f2 = f2[_R]) kinemati- 
cal fields. In fact, according to the Poincare-Wavre theorem 
(e.g., see Lebovitz 1967, Tassoul 1980), cylindrical rotation 
is equivalent to barotropicity, or to the fact that the accel- 
eration field at the r.h.s. of equations (|34[) derives from the 
effective potential 



3 ROTATIONAL FIELDS 

A natural question to ask is about the nature of the kine- 
matical fields that can be supported by the found density- 
potential pairs 5ft(\& c ) — 5ft(p c ) when considered as self- 
gravitating stellar systems. In general, the associated two- 
integrals Jeans equations (e.g., BT, Ciotti 2000) can be 
solved numerically after having fixed the relative amount 
of ordered streaming motions and velocity dispersion in the 
azimuthal direction (e.g., see Satoh 1980, Ciotti & Pellegrini 
1996). For simplicity here we restrict to the investigation of 
the isotropic case, when the resulting Jeans equations are 
formally identical to the equations describing axysimmetric, 
self-gravitating gaseous systems in permanent rotation: 

I dP - 9 ® 

p dz dz ' 

(34) 

p OR oR 

(in order to simplify the notation, in the following we intend 
p — 5ft[p c ] and $ = 5ft[<I> c ]). The quantities p, P and Q denote 
the fluid density, pressure and angular velocity, respectively; 
the rotational (i.e., streaming) velocity is v v = Q.R, while 
vr — v z = 0. 

In problems where p and 4> are assigned, the standard 
approach for the solution of equations (|34[) is to integrate 
for the pressure with boundary condition P(R, oo) = 0, 



P = 



dz ' 



(35) 



and then to obtain the rotational velocity field from the 
radial equation 



2 RdP n 3$ 



(36) 



However, v 2 , can be also obtained without previous knowl- 
edge of P, because by combining equations (|55Tl - ((3"6")l and in- 
tegrating by parts with the assumption that P — pdQ/dR = 
for z = oo it follows that 



n 2 {R')R'dR' 



(38) 



(where 7?o is an arbitrary but fixed radius), so that the gas 
density and pressure are stratified on $ c ff. 

However, barotropic equilibria are just a very special 
class of solutions of equations (|34[l . For example, in the 
context of galactic dynamics isotropic axisymmetric galaxy 
models often show streaming velocities dependent on z (for 
simple examples see, e.g., Lanzoni & Ciotti 2003, Ciotti & 
Bertin 2005). Such baroclinic configurations (i.e., fluid sys- 
tems in which P cannot be expressed as a function of p only, 
and $ c ff does not exist at all), have been studied in the 
past for problems ranging from geophysics, to the theory of 
sunspots and stellar rotation (e.g., see Rosseland 1926), and 
to the problem of modeling the decrease of rotational ve- 
locity of the extraplanar gas in disk galaxies for increassing 
z (e.g., see Barnabe et al. 2005, 2006). Note that a major 
problem posed by the construction of baroclinic solutions 
is the fact that the existence of physically acceptable solu- 
tions (i.e. configurations for which v 2 , > everywhere) is 
not guaranteed for arbitrary choices of p and However, 
the positivity of the integrand in equation (|37[1 for z > 



is a sufficient condition to have nj, > everywhere. In fact, 
in Barnabe et al. (2006) several theorems on baroclinic con- 
figurations have been proved starting from equation (J37J): 
in particular toroidal gas distributions are strongly favoured 
in order to have > in presence of a dominating disk 
gravitational field. We note that also in the present cases 
the density is of toroidal shape, altough the distribution is 
self-gravitating. 

In the case of the 5R(p c )— 5ft(\& c ) pair of the shifted Plum- 
mer sphere, the (normalized) commutator in equation (|37|l 
is a very simple function 



dp_ 9£ _ dp 9$ 
dRd? ~ Ite'dR 



15a 2 Rz(l -a 2 +r 2 
~k d7 



(39) 



which is positive for z > and < a < 1, so that v 2 -, > 
as far as 5R(p c ) is positive everywhere, i.e., for < a < a m . 
Remarkably, the integration of equation (|37[) can be easily 
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Figure 2. Isorotational contours in the meridional plane of 5R(p c ) of Plummcr sphere for a = 1/2 (top left), and a = 23/40 (top right). 
In the bottom panels contours of constant pressure for the same models. 



carried out in terms of elementary functions 
3a 2 R 2 ' 



2 



27Td 5 



6a 4 (l + R 2 



f / 4 



rf 4 (l + a 2 + r 2 ) 
12a 4 (l + i? 2 ) 3 X 



2a 2 (l + iT) 1 + a 2 + r 2 



1 + a 2 + r 2 



>, (40) 



while the direct integration of the pressure equation is 
much more complicate (again showing the relevance of equa- 
tion [37] in applications), and we report its normalized ex- 
pression on the equatorial plane only 



P{R,o) 



2(1 + R 2 



16tv(1 + R 2 ) 2 (l - a 2 + R 2 ) 2 



(41) 



Note that v 2 , oc a 2 , and so it vanishes (as a consequence 
of the imposed isotropy) when reducing to the spherically 
symmetric seed density: in fact, it can be proved that the ar- 
gument in parenthesis in equation (40) is regular for a — *• 0. 
The baroclinic nature of the equilibrium is apparent, and 
confirmed by Fig. 2, where we show the contours of constant 
vjp for a — 1/2 and a — 23/40 (top panels, correspondig to 
the models plotted in the top panels of Fig. 1). The normal- 
ized streaming velocity field in the equatorial plane is 



2 .„ n . a 2 R 2 6(l + R 
vJR,0) = — 



4o 2 (l + R 2 ) +a 4 



(l + J? 2 ) 3 (l-a 2 + J R 2 ) 3 / 2 



(42) 



and has a maximum at ii<l, while the normalized circular 
velocity is the same as that of a Plummer sphere of scale- 
lenght \/l — a 2 (cf. the comment after equation [26]), i.e., 
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l(R)=R 



d$(R,0) 
OR 



R 2 



(1 - a 2 + i? 2 ) 3 / 2 ' 



(43) 



It is easy to verify that equations (41) and (42) satisfy equa- 
tion (36). 

Unfortunately, we were unable to obtain the explicit for- 
mulae for the analogous quantities for the shifted Isochrone 
model. However, it is easy to prove that also in this case 
the commutator is proportional to a 2 , and so it vanishes (as 
expected) for a — » 0. A comparison of the rotational fields 
of the two models in the far field can be obtained by using 
the asymptotic expressions (25)-(26) and (32)-(33) can be 
used. In particular, for the real part of the shifted Plummer 
sphere we found 



2 2a 2 R 2 5a 2 R 2 (3-a 2 -9a 2 fi 2 ) ~ _„ 

v <f = — ^7 + 0{Rr ) 

while for 5R(p c ) of the Isochrone shifted sphere 



(44) 



a 2 R 2 



2 



/4 _ 13 
V3 5r 



> 



a 2 7? 2 (3159 + 5780a 2 + 21780a V 



+ 0(R 2 r- 



(45) 



Note that, even though the radial asymptotic behavior of the 
density is different in the two considered cases, the obtained 
velocity field has the same radial dependence in the far field, 
and in both cases v 2 decreases for increasing z and fixed 



R. In addition, v„ 



for r — + oo, as expected from the 



spherical symmetry of the density and potential for r — > 
oo: in fact, it is apparent from equation (|37[) that spherical 
(additive) components of p and $ (e.g., the leading terms in 
the asymptotic expansions [25]- [26] and [32]- [33]) mutually 
cancel in the commutator evaluation. 



4 CONCLUSIONS 

In this paper we have shown that the complex-shift method, 
introduced in electrodynamics by Newman, and throughly 
studied among others by Carter, Lynden-Bell and Kaiser, 
can be extended to classical gravitation to produce new and 
explicit density-potential pairs with finite deviation from 
spherical symmetry. In particular, we showed that after sep- 
aration of the real and imaginary parts of the complexified 
density-potential pair, the imaginary part of the density cor- 
responds to a system of null total mass, while the real com- 
ponent can be positive everywhere (depending on the origi- 
nal seed density distribution and the amount of the complex 
shift). 

As simple application of the proposed method we illus- 
trated the properties of the axysimmetric systems resulting 
from the shift of the Plummer and Isochrone spheres. The 
analysis revealed that the obtained densities are nowhere 
negative only for shift values in some restricted range. This 
property is linked to the presence of a fiat "core" in the 
adopted (spherically symmetric) seed distributions. It is thus 
expected that also other distributions, such as the King 
(1972) model, and the Dehenen (1993) and Tremaine et al. 
(1994) 7 = model, will lead to physically acceptable shifted 
systems. For the two new exact density-potential pairs we 
also found that the associated rotational fields (obtained un- 
der the assumption of global isotropy of the velocity disper- 
sion tensor) correspond to baroclinic gaseous configurations. 



An interesting issue raised by the present analysis is the 
toroidal shape of the obtained densities. While we were not 
able to find a general explanation of this phenomenon, some 
hints can be obtained by considering the behavior of the 
complex-shift for a — > 0. In fact, in this case it can be easily 
proved that while the odd terms of the expansion correspond 
to the imaginary part fo the shifted density (and contain the 
factor z), the even terms are real and, to the second order 



5R(p c 



2 j 2 2 

a dp a z 
2r dr 2r 3 



This expansion is functionally similar to that obtained by 
Ciotti & Bertin (2005) for the case of oblate homeoidal dis- 
tributions, where the term dependent on z 2 is immediately 
identified with a spherical harmonic (e.g., Jackson 1999) re- 
sponsible of the toroidal shape. It is clear that when restrict- 
ing to the small-shift approximation (46), all the computa- 
tions presented in this paper can be carried out explicitly 
for simple seed distributions. Moreover, from equation (46) 
it is easy to prove that negative values of density appears 
on the z axis - no matter how much the shift parameter is 
small - when the seed density has a central cusp p oc r ' 
with 7 > (e.g., see the explicit example in Appendix). 

In this paper we restricted for simplicity to spherically 
simmetric seed systems, which produce axysimmetric sys- 
tems. A natural extension of the present investigation would 
be the study of the density-potential pairs originated by 
the complex shift of seed disk distributions, such as the 
Miyamoto & Nagay (1975) and Satoh (1980) disks, or even 
the generalization of homeoidal quadrature formulae (e.g., 
see Kellogg 1953, Chandrasekhar 1969, BT). Note that in 
such cases, at variance with the spherical cases here dis- 
cussed, the shift direction is important: for example, the 
complex shift of a disk along the z axis will lead to a differ- 
ent system than an equatorial shift, that would produce a 
triaxial object. 

A second interesting line of study could be the use of the 
complex-shift to produce more elaborate density-potential 
pairs by adding models with a weight function w(a), i.e., by 
considering the behavior of the linear operator 



Qo = I p( 



i)u;(a)d 3 



(47) 



where the integration is extended over some suitably chosen 
region. We do not pursue this line of investigation, but it is 
obvious that also shifted density components with negative 
regions can be accepted in this context, as far as one is able 
to construct a weight function leading to a final distribution 
nowhere negative (as it happens, for example, in spherical 
harmonics expansions) . 
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APPENDIX A: THE SHIFTED SINGULAR 
ISOTHERMAL SPHERE 

In this case the normalized density and potential are p = 
1/r 2 and <I> = Anlnr, so that from equations ([lOj-JTT]) 



Pc 



+ i- 



2az 



(r 2 - a 2 ) 2 + 4a 2 z 2 (r 2 — a 2 ) 2 + 4a 2 z 



$ c = 27rln(r —a — 2iaz). 



(Al) 
(A2) 



Note how, at variance with the Plummer and Isochrone 
spheres, in the present case 5R(p c ) describes a density dis- 
tribution negative inside the open ball r < a and positive 
outside. The density on the surface r — a is zero except on 
the singular ring R — a in the equatorial plane. In the far 
field the density decreases as 1/r 2 , while near the center it 
flattens to — 1/a 2 , while the potential generated by K(p c ) is 
given by 



K($ c ) = 7T ln[(r 2 - a 2 ) 2 + 4a V] 



(A3) 



As for the other two cases discussed in this paper, the density 
£j(pc) vanishes on the equatorial plane, is positive above and 
negative below. 
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